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Abstract 

We give a systematic analysis of forward scattering in 3+1-dimensional quantum 
gravity, at center of mass energies comparable or larger than the Planck energy. We 
show that quantum gravitational effects in this kinematical regime are described by 
means of a topological field theory. We find that the scattering amplitudes display a 
universal behaviour similar to two dimensional string amplitudes, thereby recovering 
results obtained previously by 't Hooft. Finally, we discuss the two-particle process 
in some detail. 
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1. Introduction 



Quantum gravitational effects grow large when the length and time scales of the 
physical process of interest fall below the Planck length i v \. This well-known fact has 
the following, necessarily heuristic, explanation. If we decide to keep our coordinates 
dimensionless, the metric g^ v has dimension (length) 2 . This dimensionality can be 
made explicit by introducing some scale L, which we may take to be some length 
characteristic of the physical situation we wish to consider, and write 

g»v = L 2 g^, (1.1) 

where g^ v is dimensionless. The Einstein action 

S E [g} = J^gR (1.2) 

has the dimension of (length) 2 , that is 

S E [L 2 g]=L 2 S E [g]. (1.3) 

Now, although in the quantum theory we are supposed to integrate over all metrics, 
it seems physically reasonable to suppose that in this integral only those metrics are 
relevant whose size is comparable to that of the physical system. In other words, 
we may reasonably assume that the rescaled metric g^ u is typically of order 1. With 
this assumption, we can now indeed conclude from ( |1.3| ) that, since S E is normalized 
by Newton's constant G = i 2 ^ the dimensionless coupling constant determining the 
strength of the quantum gravitational effects is given by the ratio 

g(L) = £ f (i.4) 

of the Planck length l v i and the characteristic length L. 

Because of this reasoning, it is generally believed that to predict scattering pro- 
cesses in which the center of mass energy exceeds the Planck energy, the full theory 
of quantum gravity (or perhaps string theory) has to be invoked. However, this is 
not necessarily true, since a scattering process between two particles is characterized 
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by two energies, i. e. the center of mass energy and the momentum transfer, and even 
if the former is very large, the latter can be small. Therefore, in the physical context 
of very high energy collisions, we can regard quantum gravity as a theory with two 
coupling constants. 

This paper is an attempt to exploit this fact to start a systematic study of high 
energy scattering in quantum gravity, in the limit where the momentum transfer 
is small compared to the Planck energy. This study was motivated by recent work 
by 't Hooft, who argued that in this kinematical regime quantum gravity can be 
treated semi-classically, and furthermore displays a remarkable formal resemblance 
to string theory |L], § • Both results will be confirmed by our analysis. We will explain 
why at small momentum transfer indeed all physical modes of the gravitational field 
become classical, except for two longitudinal modes. We then show that the quantum 
dynamics of these modes is described by means a topological field theory, which 
can be characterized as a two-dimensional topological sigma model with the group 
of diffeomorphisms in 1+1 dimensions as target space. We will further study the 
interaction of this model with ultra-energetic particles and show that the amplitudes 
describing the forward scattering of these particles can be rewritten in a form similar 
to 1+1 dimensional string amplitudes. We will end with a discussion of the two- 
particle process. 



2. The Einstein Action with Two Coupling Constants 



To characterize the scattering processes that we are interested in, let us imagine 
that some day the technology will be available to build a linear accelerator that can 
accelerate elementary particles, say electrons, to energies comparable or larger than 
the Planck energy. For concreteness, let us choose our x axis along this accelerator, 
and introduce the notation x a = (t, x) for the longitudinal and y % = (y, z) for 
the transversal coordinates. Thus the momenta of the particles are of Planckian 
magnitude in the x Q -plane, while the perpendicular momenta Pi are negligible. Now, 
since the transversal cross sections of the electron beams used in our Gedanken 
experiments are necessarily very large compared to the Planck length, in practice 
all collisions between the left- and right-moving electrons have relatively very large 
impact parameters. Hence we are dealing almost exclusively with forward scattering. 
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The physics of these forward scattering processes is characterized by two length 
scales of very different order of magnitude, namely £,, , the typical longitudinal wave- 
length of the particles, and £ ± , the characteristic length of the impact parameter. 
Hence, we can form the two dimensionless ratios 

9„ = T ~ 1 

9, = £ f « 1. (2-1) 

We will now show that, in this regime, quantum gravity indeed naturally splits into 
two sectors, one strongly coupled with g , and one weakly coupled with g ± . 



2.1. Splitting the Einstein Action 

Given the above geometrical setup, it is natural to make the following gauge 
choice 

Here g a p has the Lorentzian signature, and hij is Euclidean. It is evident that locally 
the metric can always be brought in this form. This gauge may look somewhat 
unfamiliar, and indeed, it is easy to see that it is not suitable for doing perturbative 
calculations around an (asymptotically) flat background. Our interest, however, is 
in doing non-perturbative calculations in g , and for this purpose it is sufficient to 
know that the above gauge choice is non-singular around sufficiently generic metrics. 
Further, in general there may be global obstructions to choosing the metric of the 
above form, but these do not occur in our situation, because we assume that the 
topology in the x a and y l direction is trivial. 

The physical motivation for choosing the gauge ( |2.2|) is that in the kinematical 
regime of interest, the characteristic scale L. of g a p is of the order of £ p i, while the 
typical length t x of is much larger. Thus, following the line of argument given 
in the Introduction, let us introduce two dimensionless metrics g a p and hy via 

9af3 = ^ 9 a/3, 

hij = f h ir (2.3) 
We again will make the plausible assumption that g and h are of order 1. 
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Now, a straightforward calculation, outlined in the Appendix, shows that the 
Einstein action in the gauge ( |2.2| ) splits into two terms 



S E \g] = S // \g,h] + S 1 .[h,g), (2.4) 

given by 



and 



S„ [g, h}= Jy/g (VhRh + \Vhh^d iga pd ]gi6 e^e ps ) , (2.5) 
S±[h, g] = JVh (yfgRg + l^gg^d^dph^ 1 ) . (2.6) 



The crucial property separating both terms is that they behave differently under 
the constant rescaling (|2.3|) of g a p and hij, namelyQ 



S^gJlh] = i'S^g, h), 

S x [£lh,£lg] = £ 2 ± S ± [h,g}. (2.7) 

Thus, after we make Se dimensionless by dividing by Newton's constant G = £ 2 vh 
we see that the S part of the Einstein action is strongly coupled, with coupling 
<?// = £pi/£/n whereas the coupling constant g ± = £ p i/£ ± of the other term S± is very 
small. We arrive therefore at the important conclusion that as far as S± is concerned 
we are in the classical regime. In the following we will exploit this simplification of 
the dynamics of quantum gravity to analyze the theory to leading order in g ± , while 
keeping track of the full non-perturbative dependence on g... 

To complete the picture we should include the ghost action associated with our 
gauge choice g^ = 0. Under an infinitesimal coordinate transformation we have 
figia _ fovQ^a _|_ g al3 dp!; 1 , so we find that the ghost action consists of the two terms 

Sgh [b, c] = S f/ [b, c] + S± [b, c] 

SAb,c] = f^gVhh^b ia d jC a , (2.8) 



S± [b,c) = JVh^gg aP b ia dpc\ 



Again we observe that after the rescaling (|2.3| ) the first term S., is strongly coupled 
and the second S ± weakly. 

*The e-symbols in ( |2.5[ ) and ( |2.6| ) contain a factor of (y^ff) -1 resp. (^/h)^ 1 , so that they 
transform as tensors. 
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2.2. Isolating the Relevant Gravitational Modes 



It is instructive to make an inventory of the different degrees of freedom of the 
gravitational field in relation with the above decomposition of the action. As usual in 
quantum gravity we have, after gauge-fixing, 6 remaining components of the metric 
and 4 anti-commuting ghost variables. The ghosts eliminate the 4 longitudinal 
modes of the metric via the BRST-invariance of the action and effectively impose 
the 'Gauss law' constraints 

R ia = 0. (2.9) 

This leaves 2 physical degrees of freedom constituting the gravitational radiation 
modes. A useful parametrization of the metric variables, which manifestly exhibits 
the different modes, is 

g a( s = e^ x ^d a X a df,X b r] ab (2.10) 

hij = exWdiY'djY'Sn. 

Here X a (x,y) and Y p (x,y) are the diffeomorphisms that relate g a p resp. to the 
diagonal metrics. As such, they are paired with the ghost fields c a resp. c l via the 
BRST symmetry. Hence they are the longitudinal modes and locally decouple from 
physical processes; as we will see, however, their global variations are important 
observables. The conformal factors and x are physical, BRST-invariant fields and 
represent the two transversal components of the graviton. 

To see which of these modes are physically relevant in our context, let us consider 
the coupling to matter. In general, we can include matter by adding a source term 
to the Einstein action 

S E [g]^S E [g] + S m [g}. (2.11) 

This source term represents the stress-energy of the matter particles, of which the 
transversal components are assumed to be negligible compared to longitudinal 
components 

1 5S m [g] 

Va $g a/3 

Here T a p is effectively traceless, provided the rest-masses of the matter particles 
are very small compared to the Planck mass. The remaining two components of 
T a/ 3 can be identified with the momenta of the particles, which are either left- or 
right-moving. We see therefore that the matter stress energy only couples to the 
off-diagonal variations of g a/ 3, which are parametrized by the fields X a . 



Tap = (2-12) 
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An important property of S ± is that it is invariant under reparametrizations x a — > 
x a (x, y). As a consequence, one finds that upon inserting the parametrization ( |2.10 ) 
into ( |2.5|) and (|2.6| ) that S ± does not depend on X a . Now, since g ± <Kl, the functional 
integral over metrics is strongly peaked near the space of field configurations that 
minimize S±, and thus to leading order in g ± , we can restrict the configuration 
space of gravity to this vacuum subspace. Combined with the above observation 
that the matter only couples to the X a fields, this shows that all other variables 
can be treated classically. The X a are therefore the relevant quantum mechanical 
modes, since their dynamics is dictated purely by the strongly coupled S,, piece of 
the action. 



3. The Topological Effective Theory. 



The restriction to the vacuum subspace of S x eliminates the local gravitational 
degrees of freedom: only the global variations of the remaining longitudinal modes 
X a can be physically relevant. In this section we will construct the topological field 
theory that describes their dynamics. 

From the above discussion we learn that, to leading order in g ± , we can restrict 
ourselves to the special class of vacuum field configurations such that 

S ± [h,g) = 0. (3.1) 

This leads to the equations 

d a hij = 

R g = 0, (3.2) 

which single out the vacuum solutions in which no graviton waves propagate. These 
solutions are of the general form 

Kj = hij{y) 

Q a p = Vabd a X a d p X b (3.3) 
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The variables X a are maps of the two-dimensional x-plane onto itself and vary in the 
transversal y directions. Here the transversal metric hij(y) is allowed to be arbitrary 
and will be treated as a fixed classical background. 

Along with the metric, the ghost fields are also constrained to solve the equation 
of motion of the weakly coupled S± piece of the ghost action ( j2.8| ). This equation 
for the anti-ghost, V a b ia = 0, is solved by 

b ia = tapd^bi, (3.4) 

where the new anti-ghosts 6j are unrestricted field variables. Substituting the 
parametrizations ( |3.3|) and ( |3.4|) into the Einstein action gives the effective action 
for the new variables X a , c a and b l . 

To write this reduced action, it is convenient to introduce the vector fields VJ" 
via the relation 

d t X a + Vi a d a X a = 0. (3.5) 

These variables Vj 01 describe the flow of the X a fields in the ^/-direction, and are 
analogous to the fluid velocity in fluid mechanics. The transversal derivatives of the 
metric can be expressed in terms of as 

d ig a ? = V a Vi P + V p Vf. (3.6) 

Here and in the following we keep the X a dependence implicit in g a p and ' . After 
substituting this into ((2.5|) and using the identity e^e^ = e ai eps — e^e^s we find 
that the action naturally splits in two contributions 

S = 5*i + 5*2 

51 = JVhV9(Rh + e aj epsV a VfVW tS ) (3.7) 

5 2 = JVhVd^-W^V/f + Veapff'diC?). 

Notice that, when we insert ( |3.3|) and (|3.5| ), this action is highly non-linear in the 
fundamental fields X a . It describes, as promised, a topological field theory: there 
are no local physical degrees of freedom. The first term is invariant under local vari- 
ations of the A a -fields, and can in fact be written as a total derivative. The second 
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term in (|3.7|) , when regarded as a two-dimensional field theory in the transversal 
y-plane, represents a kind of two-dimensional topological sigma model [|]]. In this 
interpretation the target space of the model is given by the space of diffeomorphisms 
of the x a plane. 

Both terms in the action ( |3.7| ) are invariant under the nilpotent BRST symmetry 

5X a = ec p dpX a , 
5V = eie^V^ + c /3 d f3 b i ), (3.8) 
Sc a = ec p d p c a , 

under which Vi° transforms as 

5V? = e(-d iC a + c?d f) Vi a - Vi p d p c a ). (3.9) 

This nilpotent symmetry removes all local fluctuations of the X a -field, and is a 
remainder of the BRST-invariance of the original theory, as can be seen as follows. 

The original anti-ghost field b{ a transforms under the BRST-symmetry into the 
Ri a components of the Ricci-tensor, which are the constraints associated with our 
gauge choice. Therefore, the redefinition (|3.4j) of the ghost field is compatible with 
BRST-invariance only if R ia is also of the form 

Ria = e^dPRi. (3.10) 

Using that the background satisfies the field equations of S±, which gives that 



V a Ri a = —V J Gij = , one deduces that (|3.10D indeed holds. Explicitly, one 
has 

Ri = t a pd a Vf. (3.11) 



The BRST-symmetry of the reduced action ( |3.7|) thus directly follows from that of 
the gauge fixed Einstein action. It effectively imposes the constraint Ri = 0, which 
tells us that the flows V{ are curl free. 

Let us end this subsection with one final comment. The attentive reader may 
have worried a little about the fact that the changes of variables from g a p to X a and 
from b ia to b l , both require some Jacobians to correctly treat the functional measure. 
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It turns out, however, that these two Jacobians cancel each other. Explicitly, if we 
define the measures of g a p and X a through 

\\6X a \\ 2 = JVhV9VabSX a 5X b , (3.12) 
one obtains from ( 2.1(J| ) that the measures of g a p and X a are related via 



[dg Q p] = [dX a ]\det(e a ^)\ 2 . (3.13) 
Similarly, since the ghosts are anti-commuting, we deduce from (|3.4j ) that 

[dbia] = [rf6 l ]|det(e Q/3 V /3 )|- 2 . (3.14) 
So the total Jacobian of the substitution (g a /3, bi a ) — > (X a , b l ) is indeed equal to 1. 



4. The Action for the Boundary Values 

Let us now ask ourselves what physical degrees of freedom are described by the 
action (|3.7| ). For definiteness, we will from now on consider the situation where the 
^-coordinates parametrize (a part of) the two-dimensional Minkowski plane M. 
Thus our space has a boundary DM in this direction, consisting of four asymptotic 
regions, two belonging to the past null infinity Z~ and two constituting the future 
null infinity X + (see fig. 1). Now we notice that local variations of the fields X a 
are not physical, since they correspond to diffeomorphisms and are cancelled by the 
BRST-symmetry. Therefore the only possible physical degrees of freedom are the 
boundary values 

X a = X a . (4.1) 

\dM y ' 

We will allow these boundary values X a to be arbitrary, because, as we will discuss 
below, this enables us to incorporate matter into the theory. Furthermore, to make 
them into true physical observables, we impose the usual restriction that gauge 
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Fig. 1: The boundary dM of the two-dimensional parameter space M of the 
longitudinal coordinates x a has four components corresponding to the four 
asymptotic regions of the Minkowski plane. 

transformations of the theory, i. e. diffeomorphisms, asymptotically go to the identity 
transformation. After gauge fixing this translates into the asymptotic condition that 
the c-ghost vanishes on dM, and thus the X a are BRST-invariant. 

We will now show that for physical field configurations the action (|3.7| ) can be 
expressed purely in terms of X° ' . This is most evident for the first term Si[X]: it is 
easy to see that it is the integral of a total derivative. Hence 

S X [X\ = S 9M \X\. (4.2) 

One finds 

S dM [X] = fjx a JVhe ab (R h X a dJ( b + dpC a d a d l X b ). (4.3) 

The second term of the action (|3.7| ) plays the role of the gauge fixing term to the 
action, needed to eliminate the redundancy caused by the topological nature of the 
first term. Following a standard argument, often applied in topological field theory 
||, we can use the symmetry ( |3.8|) of the action to restrict the functional integral 
to the BRST-invariant field configurations. This implies that we may put the ghost 
fields c a equal to zero, and in addition impose the constraint 

e a pd a Vf = ^d a {diX a dpX a ) = 0. (4.4) 

When we insert this into fl3.7| ) we find that the SVpiece of the action vanishes 
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identically. We conclude, therefore, that for physical field configurations the total 
action indeed reduces to the simple boundary term (|4.3|) . 

Finally, a more precise specification of what we mean by the asymptotic values 
X a of X a is perhaps in order here, since a priori at least one of them becomes infinite 
when we approach the null infinities X^. The point here is that, as we will see, only 
variations of the asymptotic values of the X a are physically relevant. Hence it is 
sufficient for our purposes to consider, instead of the full Minkowski plane, only 
a finite region M, large enough so that variations 5X a at the boundary dM are 
equal to the asymptotic variations. This procedure will be implicitly adopted in the 
following sections. 

4.1. Coupling to Matter 

We will now discuss the coupling of external matter to the gravitational field, 
and show that it can also be represented entirely in terms of the boundary values 
X a . We consider the situation in which asymptotically all particle momenta are in 
the ^"-direction, so that at the boundary dM all components except T a p vanish. 

First let us assume that also in the interior region T a p are the only non- vanishing 
components. Its conservation law then reads 

V Q T Q/3 = 0. (4.5) 

This tells us that the matter source term S m [g] is invariant under diffeomorphisms 
of x a . In terms of the X a field variables this implies that S m [g] is invariant under 
local variations, and indeed only depends on the boundary values X . To write 
this dependence, it is convenient to represent the incoming and outgoing matter 
at the boundary dM by means of a momentum flux P aa , defined in terms of the 
stress-energy tensor via 

T aP = P aa dpX a . (4.6) 
Then the variation of the matter contribution with respect to X a is 

8S m [g(X)\ = 2 JVhV9V a (P a JX a ) 

= 2<ldx a fy/hej P a0 5X a . (4.7) 

J dM J 
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Thus, the matter momentum in- and out-flux couples directly to the variations of 
the boundary values X a . 

In general, however, it is possible that in the interior of M also the mixed compo- 
nents Ti a of the stress tensor become non-vanishing. Although in this case we have 
to refine the above derivation of the matter coupling, it in fact turns out that the 
final result ((O]) remains unaltered. Let us indicate how this goes. The conservation 
equation ( |4.5| ) is now modified to 

V a T al3 + V% = 0. (4.8) 

Because we assume = 0, we know that T ia in addition satisfies V Q T ia = 0. It is 
therefore of the form 

T ia = eap&Ji. (4.9) 

Here J, can be identified with the current associated with the boosts SX a = e%X b in 
the X-plane. Namely, by integrating ( |4.8|) over the plane M against d^X a X b e a b 
and inserting the redefinitions (|4.6|) and Q4.9|) we deduce that for all y 

V, f Jgr = I dx a eJPSX b e ab , (4.10) 

JM JdM ' 

which expresses the conservation of angular momentum. 

When Ji is non- vanishing, the matter action S m [g] is no longer invariant under 
local variations of X a . However, coordinate invariance tells us that this local X a - 
dependence must cancel when we include the gravitational interaction. To see this, 
consider the BRST- variation of the matter action S^fg]. It can be expressed in 
terms of the current Ji as follows 

SbrstSmig] = 2e j ' Vh^/gT a pV a cP 

= 2e jVhy/gfe^dPdiC^ (4.11) 

where we used ( |4.8| ) and (f4.9|). Thus the combined matter and gravitational action 
is BRST-invariant, provided we include J 1 in the transformation law (|3.8|) of the 
anti-ghost. As a consequence the constraint fl4.4j) on physical field configurations is 
replaced by 

e af sd a V^ = Ji. (4.12) 
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With this modification we can now follow the argument of the preceding subsection 
to show that, under the above physical constraint, the dependence of S m [g] on local 
variations of X a indeed cancels against that of the second term of the gravitational 
action (|3.7|). Hence also in this case the entire coupling to matter is described by 
equation (|4.7p . 



5. The Scattering Matrix 

In this section we will determine the amplitudes describing the gravitational 
scattering of the matter particles. First we like to summarize the result of the 
above analysis. We have shown that, due to the BRST-invariance, the functional 
integral over the fields X a , b l , and c a is concentrated on the space of physical field 
configurations, satisfying the constraint ( |4.12j ). For given boundary values X a , this 
constraint uniquely fixes the interior values of X a , and thus all dynamics is reduced 
to the boundary DM. It will be convenient to introduce a 'time'-variable r which 
parametrizes the coordinates x a (r) on DM. Notice that, since DM is closed, r is 
periodic. The action that describes the dynamics of the boundary values and its 
coupling to matter is 

S dU [X] = S[X] + 2jdrjVh P a , T X a , (5.1) 

with 

f f BX a 

S[X] = JdrJ^fh e ab — (A h - R h )X b , (5.2) 

where A/, denotes the scalar Laplacian in the transversal ?/-plane. Here we used 
(|4.7|) to represent the matter coupling in terms of the momentum density P a r = 
epP aa d T xP . In (O) and from now on we drop the bar on the boundary values of X a . 



In the rest of this section it is understood that all the fields are defined on DM. 



5.1. Quantum mechanics of the boundary values. 

The boundary theory (|5.1| ) can naturally be regarded as a (non-relativistic) 2+1- 
dimensional field theory. Indeed to obtain the scattering amplitudes we will still 
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need to perform the remaining functional integral over the X a , and it is therefore a 
natural step to consider the quantum mechanics of this boundary theory. For this it 



is a fortunate fact that the action (|5.2| ) is quadratic in the fields X a } especially since 
it was derived from the highly non-linear Einstein action. Canonical quantization is 
therefore straightforward. It leads to the following commutation relations 

[X a (y l ),X b (y 2 )]=ie ab f(y 1 ,y 2 ), (5.3) 

where 

(A h -R h )f(y 1 ,y 2 ) = 6^\y 1 ,y 2 ) (5.4) 

defines the Green function of the modified Laplacian (A^ — Rh). This commutator 
previously appeared in the work of 't Hooft 0. Intuitively it can be understood from 
the Einstein equation (|5.9| ), which relates the coordinate fields X a to their canonical 
conjugate momenta. Let us denote by 7i the Hilbert space in which the commutator 
algebra ( |5.3|) is realized; vectors in 7i can be represented as functionals of X + (y) 
or of X~(y). Note that, since the action is invariant under reparametrizations in r, 
there is no Hamiltonian, i.e. H = 0. 

The matter coupling is described in terms of 'vertex operators' of the form 

V(P) = Pexp(2zJdrJ^P^ T (T,y)X a (y)). (5.5) 

Here the path-ordering is necessary, as the X a (y) do not mutually commute. The 
momentum flux P a ^ T will for the moment still be treated as a classical quantity, 
but in the next subsection, when we study the scattering matrix, P a ^ T becomes a 
quantum mechanical operator in the Hilbert space of the matter theory. The idea 
is that in this way, by first computing the expectation value of the general vertex 
operator (|5.5p, we obtain a generating functional for all scattering amplitudes. Since 
the boundary time r is periodic, this expectation value is given by the trace over 
the Hilbert space Ti 



V{P)) = tT n {Pexp[2i^drJVhP atT X a ]). (5.6) 



The result can be computed from the commutation relations ( |5.3|) via a simple 
application of the CBH-formula, or alternatively, by using the relation with the 
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functional integral. One finds 



V(P)) = ^[ijdrdr' jVhe ab P a>T (A h - R^P^]. (5.7) 

In addition we obtain a delta-function imposing the condition that the total mo- 
mentum flux at the transversal position y vanishes 

<fdrP a ^r,y) = 0. (5.8) 

This condition, which is due to the invariance of S[X] under r- independent shifts 
of X a , guarantees that the result ( |5.7| ) does not depend on the base-point of the 
r-integration. 

The result ( |5.7|) for the generating function is exactly given by the semi-classical 
expression, exp(iS[X c i]), where X c i solves the classical equation of motion 

f)X a 

(A h -R h )— = e ab P b , T . (5.9) 

This equation, which is in fact the remainder of the Einstein equation, only involves 
the values of the coordinate fields X a and momenta P a on the boundary dM. Thus 
with our formalism we have succeeded in computing the Einstein action of the 
relevant classical solution, without needing to integrate the equation of motion for 
the interior region. 

In principle the classical configuration for X a in the interior region is determined 
by the physical constraint ( [4.12 ) involving the current Jj. The integral of ( 4.12 ) over 



the two-dimensional parameter space M gives for all y the relation 

r f)X a r 

f dr— d.Xaiy) = / y/gJ^y). (5.10) 
J 8M or Jm 

where we used that e^d 01 Vf = e a Pd a (diX a dfsX a ). As a consistency check one can 
verify, using the equation of motion ( |5.9| ), that this relation implies the conservation 
law gTTCp . 
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Fig. 2: The amplitude ( |5.13| ) is the S-matrix element between the state with 
momentum flux P™ at the in-region X~ and the state of given out-going flux 
P° ui at J+. 



5.2. Description of the Scattering Matrix 

Let us now determine the scattering matrix for (point-like) matter particles. 
Since we assume their masses to be negligible compared to the Planck mass, they 
are divided into left-movers with momenta p l _ and right-movers with momenta p 3 + . 
Now, it is important to note that the left- and right-movers come in from or go 
out to different asymptotic regions of space-time, so that the boundary DM is di- 
vided into four different components, on each of which either P_ T or P +)T vanishes, 
(see fig. 2). Correspondingly, the vertex operator (|5.5| ) representing the momentum 
flux decomposes into the product of four vertex operators, where, for example, the 
operator representing the right-moving incoming particles has the form 

V(P?) = exp(2* J VhPl n (y)X + (y)), (5.11) 

where 

P l + n (y) = JjrP + , T (T,y) (5.12) 

is the total incoming momentum at the transversal position y. The scattering am- 
plitude is obtained by specializing the general formula ( |5.7|) to the present situation. 
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This gives the result 



V(P; n ) V{P™) V*(Pl ut )V*(P mt )) = S(P) 5(P° a ut - Pf), (5.13) 
where 

S(P) = exp(z JVhJVhr ] ab P a (y 1 )f(y 1 ,y 2 )P b (y 2 )) (5.14) 

is the scattering matrix. 

The variables P a {y) act as quantum mechanical operators in the matter Hilbert 
space, and can be used to characterize the in- and out states. Because the operators 
P a {y) mutually commute, we can find a basis of eigenstates in the in- and out- 
Hilbert space. Such a basis is obtained by specifying for each particle its longitudinal 
momentum and transversal coordinate y^K On the iV-particle state \p^\y^), 
the momentum operator P a {y) has the eigenvalue 

N 

P a (y) \pf,y {i) ) = £^(y,y (i) ) \ P (l) ,y {l) ). (5.15) 

i=l 



Since the S- matrix ( |5.14j) is expressed purely in terms of the P a {y), it is also diagonal 



in the same basis, with eigenvalues 

S(p®,y®) =J[exp(iT } *p®pPf(y®,yV))). (5.16) 

Equations ( |5.14j) and (|5.16|) give the final result for the scattering matrix. Notice 
that the iV-particle amplitude factorizes into the product of two-particle amplitudes, 
which tells us that there are no three- or higher point interactions. In the next section 
we will discuss the two-particle process in more detail. 

We observe that the 5-matrix element ( 5.16Q can be written in a form similar to 
a string amplitude as follows 

Stp®,y®)=[[dX]e iS '»W JJ e wrf^»(« w ), (5.17) 

i 

where S str denotes the two-dimensional action for a string embedded in a 1 + 1 
dimensional space-time (cf. 0) 

S str [X] = J (PyVh(^ j d i X a d j X a + R h X a X a ) . (5.18) 
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An unusual feature is the factor i in front of the action, even though the 'world 
sheet', the y-plane, has Euclidean signature. 

Another difference with the usual string amplitudes is that there is no integration 
over the positions of the vertex operators, nor over the metric hy. The reason for 
not integrating over these variables is that we assume that for large x a we are in 
the classical regime, where general covariance is broken in the ^/-direction via some 
vacuum expectation value (hij) = h^. One can imagine, however, some physical 
situations where general covariance remains unbroken. In this case we do have to 
integrate over and the formal analogy with string theory becomes almost exact. 



6. The Two-Particle Amplitude 

In this section we will take a more detailed look at the two-particle scattering 
process, and discuss its geometrical interpretation. In the following we will assume 
that the transversal metric is flat = 5ij. In this case the Green function / (2/1,2/2) 
becomes 

2/2) = — log — 2/2 1 2 - (6.1) 



6.1. The Shock Wave Geometry 

The scattering matrix obtained in the previous subsection has a concrete physical 
interpretation in terms of gravitational shock waves H . To exhibit this interpre- 
tation, we first observe that the N-particle amplitude factorizes into the product of 
two-particle amplitudes, given by 

S12 = exp(2ipip£f(y 1 ,y 2 )) (6.2) 

where pf and 2/i are the momentum and transversal coordinate of particle i. This 
two-particle result is obtained from ( 5.14Q by choosing the in-flux of momentum at 



X of the following form 

p irX x ^y) = Pi 8( x+ - x i) 5 (y-yi)i 

p in( x iv) = pt 5 ( x ~- x 2) s (y-y2)- (6-3) 
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x =X 



Fig. 3: The shock wave geometry of a fast left-moving particle with momentum 
Pi and transversal position y\ has for all y a discontinuity at the ^"-trajectory 
of the particle. 



The corresponding classical solution of (oTa) is given by 

x in = x ~ ~ PiO(x + -x+) f(y,yi), 

Xfn = % + - Pt6(x--x 2 ) f(y,y 2 ). (6.4) 

These classical field configurations for Xf n and X~ n both exhibit a discontinuity at 
the x a trajectory of one of the particles. These discontinuities are the well-known 
gravitational impulsive wave, or 'shock wave', solutions to the Einstein equation 
obtained in The form of this shock wave geometry for a single particle is in- 
dicated in fig. 3. It shows that a right-moving geodesic is shifted by an amount 
Pif(y,yi) when it passes through the x Q -trajectory of the fast left-moving particle 
with momentum p\. 

The formula ( |6.4j ) only describes the two gravitational shock waves in the in- 
region, i.e. before they collided. It is clear that in the out-region the solution will 
become more complicated. To solve for the complete metric describing the two col- 
liding shock waves indeed appears to be a rather non-trivial problem. The technical 
difficulty here is that after the two particle trajectories cross each other, the posi- 
tion of each shock wave is shifted by the discontinuity of the other. This distorts 
the waves in such a way that they are no longer stable solutions of the Einstein 
equation. What we know, however, is that the out-flux P Ta (r,y) at X + must, by 
momentum conservation, yield the same integrated momentum distribution P a {y) 
as the in-flux. This is all information we need to compute the scattering amplitude 
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S12 via the general formula (|5.14 ). 

Finally, let us note that the amplitude (|6.2| ) has been derived in an alternative 
way in by considering the effect of the shock wave of one of the two colliding 
particles on the wave packet of the other. This reasoning makes use of the fact that 
one can go to a frame where one of the particles moves slowly, so that only the 
gravitational field of the other fast moving particle needs to be considered. As an 
operator in the Hilbert space of the slow particle, the S'-matrix ( |6.2| ) can then be 
recognized as the generator of the coordinate shift caused by the shock wave of the 
fast particle. 



6.2. Partial Wave Analysis 

The two-particle process is special because in this case our geometrical setup in 
fact does not break Lorentz invariance: the division into longitudinal and transversal 
directions puts no restrictions on the momenta of the two particles, but is defined 
in terms of them. We will now indeed show that the two-particle S-matrix and 
differential cross section can be given in a Lorentz invariant form. Our approach is 
based on a partial wave analysis of the scattering amplitude, and can be regarded as 
a refinement of the method of QIJ . It was found in |] that the two-particle amplitude 
possesses an infinite number of poles for imaginary P 2 , which were conjectured to 
indicate the presence of bound states. One of our aims is to determine whether 
these poles already occur at the level of partial waves or not. 

To write the S'-matrix ( 6.2| ) in manifestly Lorentz invariant form, we can first 
replace P1P2 by P\P2p,- Next, let us introduce the total linear and angular momentum 
operators 

P" = P l L+p£, (6-5) 

JfJ. = |p| e /i^Ao-PlP2 (^1 — X V)i 

and express PiP2 l i and the impact parameter \y\ — 2/ 2 1 in terms of these quantities 
as (here we use that p\ = p\ = 0) 

2rfp2 M = P\ (6.6) 
\yi-V2? = J 2 IP 2 - (6.7) 
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Notice that this implies that as a quantum mechanical operator the impact param- 
eter has a discrete spectrum. 

The two-particle scattering operator is conveniently expressed in terms of the 
mutually commuting operators P and J as follows 

(J 2 \~ ip2 

Si2 = [p- 2 J ■ (6.8) 

Hence the S-matrix is diagonal in the partial waves basis of eigenmodes \P,l,m) of 
P, J 2 and J 3 , and can be represented as 

S l2 \P,l,m) = e 2i5l{p) \P,l,m), (6.9) 

where the phase shift Si(P) is given by 

6 l (P) = - 2 -P 2 \o g (^±^). (6.10) 

This result gives the exact leading behaviour of the scattering matrix for large 
angular momentum I, that is 

/(/ + 1)>P 2 , (6.11) 

with P 2 measured in Planck units. Notice therefore that if P 2 is small our approx- 
imation is reliable for all partial waves, except for the s-wave. Indeed, the above 
result for the phase shift becomes singular for 1 = 0. Notice further that, for fixed 
I, the S'-matrix does not exhibit any poles in P 2 . 

From the knowledge of the phase shifts Si(P) one can compute the differential 
cross-section via 

> = \m? 

1 oo 

fptf) = ^£(2/ + l)(e^-l)fl(cos0). (6.12) 



Inserting the result ( |6.10| ) (and restricting the sum in ( |6.12| ) to I > 0) we obtain an 



expression for the amplitude fp{0) which is reliable for small scattering angles 9. 
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To extract this leading behaviour we make use of the following asymptotic formula, 
which becomes exact for large I and small 9 



Pi(cos6) -»• J o ((2/+l)sin|0), (6.13) 
where Jq is the first Bessel function. Let us now re-introduce the impact parameter 



y = y 1(1+1) /\P\ and approximate the sum over I by the integral over y. The 
resulting integral can be performed and gives in this approximation the following 
leading behaviour 

f p (Q) = -i\p\ J^dyy 1 - 21 ? 2 J (2y\P\ sin \0) (0 small) 

2sin 2 ^r(l + *P 2 ) ^ S 2 ^ • (b ' i4j 
This formula agrees with the result obtained in |TJ. 



The most striking feature of the scattering amplitude ( p.!4[ ) is the occurrence 



of poles at the center of mass energies P 2 = —iN, in Planck units. However, it is 
clear from the above derivation that they do not correspond to any resonant states 
of quantum gravity, but are produced by the replacement of the sum over I by the 
integral over y. If there would exist resonances with definite spin, the corresponding 
poles would already have been visible in the phase shifts Si(P). The poles in ( |6.14| ) 
come from the y — > part of the integral and disappear when this part of the integral 
is cut-off and replaced by a sum. Hence, assuming that the s-wave contribution has 
no poles, the original expression Q6.12 ) for the scattering amplitude is analytic for 



all values of P 2 . We conclude therefore that result ( |6.14| ) describes the correct small 
9 behaviour of the amplitude as long as ImP 2 > —1; beyond this point subleading 
terms in ( |6.12[ ) become important. 



7. Concluding Remarks 



In this paper we have described a systematic method for studying forward scat- 
tering in quantum gravity. It exploits the presence of the small dimensionless pa- 
rameter g± = (l p i/l ± ) 2 , in terms of which one can do perturbation theory. The final 
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result Q5.14 ) for the scattering amplitude is valid only to leading order in g±, but is 
non-perturbative in the other dimensionless coupling constant g.. = (lpi/1^) 2 . It can 
therefore be applied for arbitrarily large center of mass energies. 

We have shown that the relevant gravitational modes that mediate between the 
high energy particles are described by means of a topological field theory and that 
their dynamics and the coupling to matter is formulated in terms of the boundary 
values of X a . This result is universal in the sense that it only depends on general co- 
variance and dimensional analysis, and also applies in the presence of a cosmological 
constant or curvature squared terms. Curvature squared terms will, however, mod- 
ify the boundary action (|4.3| ) by higher derivative terms. The fundamental property 
that remains is the absence of local dynamics at short distances. This is a reflection 
of the fact that in the ^"-direction general covariance in unbroken, and in this sense 
we have truly dealt with the 'topological phase' of quantum gravity. 

The fact that the leading order amplitude (|5.14 ) is unitary, tells us that to this 
order no gravitational radiation is produced in the scattering process. Indeed, the 
coupling to gravitons is of lower order in the momentum transfer. It is possible to 
compute these higher order corrections due to the presence of gravitons, although 
at some point one is bound to encounter the problem of non-renormalizability of 
quantum gravity. String models of quantum gravity are likely to avoid this problem, 
but we expect that these string modifications will not drastically alter the leading 
order result (see, however, ||). 

One of the more intriguing conclusions of our paper is the commutation relation 
( |5.3| ) between the coordinate fields X a . Its geometrical origin lies in the fact that the 
discontinuities describing the gravitational shock waves ( |6.4|) of the left- and right- 
moving particles are in a way incompatible: the shift caused by the left-mover causes 
an uncertainty in the location of the right-mover, and vice versa, at the point where 
the two x a trajectories intersect. One may speculate that the commutator hints 
towards the existence of a new quantum gravitational uncertainty principle, which 
could have radical consequences. Indeed, 't Hooft has proposed that, by taking this 
uncertainty principle serious, one can construct an S'-matrix for the successive black 
hole formation and evaporation process 0. 

In view of this application, it would be very interesting to generalize our analysis 
to the case where x, instead of being a cartesian coordinate, represents a radial 
coordinate r. This corresponds to the physical situation in which the momenta of 
all particles are directed to the same point (r = 0) in space. However, here we 
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will meet a new difficulty, since for very small values of r our assumption about 
the transversal momenta being small becomes unattainable. Therefore additional 
modes of the gravitational field need to be included in the analysis. One indeed 
expects that in this situation black hole formation may take place, and then e.g. the 
mass of the black hole becomes a dynamical variable. 



Appendix: The Einstein Action in the Gauge (2.2) 



In this Appendix we outline the derivation of the formula ( |2.5| )-( |2l]| ) for the 
Einstein action in the gauge 

*- - ( 7 I, ) ^ 

Here and in the following fi, v run from to 3, a, /3 from to 1 and i,j from 2 to 3. 
Under this decomposition the Ricci scalar of g^ v reads 

R W = gOfigTR^ps + Wh kl R tkjl + 2g a ?h i m am (A.2) 

where R^ u \ a is the Riemann tensor 

Inserting the explicit form of the following Christoffel symbols 

^ aij ^ija ~^f^ot\^ij (-^-'^O 

while keeping r^g and rf- implicit, we find for the first two terms in (|A.2|) 
g^g^R a ^8 = Rg- \h ij d iQaP d j9lS (g ai ; g^ - g<# <g^ s ) 
h i3 h kl R ikjl = R h - \ g a/3 dahij dph kl (h a V k -h l ih kl ) (A.5) 
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where R g and Rh are the scalar curvatures of g and h considered as two-dimensional 
metrics. The last term in ( |A.2 ) becomes 



2g a P0R odPj = | h ij d i9aP d j9lS (g aS g^-g af3 g^) 

+ \ g af3 d a hij dph k i {h il h jk -h ij h kl ) (A.6) 

Vg Vh 

Here the last two terms contribute total derivatives to the Einstein action, and can 
therefore be omitted. Adding the remaining terms and using 

h u h jk_ h ij h ki = e ik e ji (A.7) 
leads to the result ( p.5|) -( p^6|) for the Einstein action. 
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